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Abstract 

We discuss the analytical solution of the two-loop sunrise graph with arbitrary non-zero 
masses in two space-time dimensions. The analytical result is obtained by solving a second- 
order differential equation. The solution involves elliptic integrals and in particular the so- 
lutions of the corresponding homogeneous differential equation are given by periods of an 
elliptic curve. 



1 Introduction 



The two-loop sunrise graph with non-zero masses is the simplest Feynman integral, which cannot 
be expressed in terms of multiple poly logarithms. It has already received considerable attention 
in the literature [fTI-fTOl. The analytical solution for the equal mass case has been discussed in |I6|. 
Up to now, less is known for the unequal mass case. The state of the art in the unequal mass case 
can be summarised as follows: First of all, expansions around special points, like zero momentum 
squared, threshold or pseudo-thresholds are known [|11-h16] . Furthermore, representations of the 
full integral in terms of Lauricella functions [2-44iT7l or one-dimensional integral representations 
involving Bessel functions [|7l[8l are also known. For practical purposes, numerical evaluations 
are available [[T8ti20l. 

The two-loop sunrise integral with non-zero masses is relevant for precision calculations 
in electro-weak physics [3J, where non-zero masses naturally occur. In addition, the two-loop 
sunrise integral with non-zero masses appears as a sub-topology in many advanced higher-order 
calculations, like the two-loop corrections to top-pair production or in the computation of higher- 
point functions in massless theories [121 [1221 . It is therefore desirable to understand this integral in 
detail, as such an understanding will pave the way for more complicated processes. As already 
mentioned, the two-loop sunrise integral provides the simplest example of an integral which 
cannot be expressed in terms of multiple polylogarithms. Therefore this integral is the ideal 
playground to learn more about transcendental functions occuring in Feynman integrals beyond 
multiple polylogarithms. 

In this paper we discuss the two-loop sunrise integral with unequal non-zero masses in two 
space-time dimensions. The restriction to two space-time dimensions allows us to focus on the 
essential difficulties and avoids some entanglements with complications which we know how 
to handle. In two space-time dimensions the integral is finite and depends - as we will see 
later - only on the second graph polynomial, but not on the first graph polynomial. Working 
in two space-time dimensions removes therefore the issues of ultraviolet divergences and the 
dependence on the first graph polynomial. Of course we are interested in the end in the result 
around four space-time dimensions. Using dimensional recurrence relations [|23ll24l . the result 
in two space-time dimensions can be related to the result in four space-time dimensions. 

Recently, it emerged that methods from algebraic geometry can be very useful for the com- 
putation of Feynman integrals. For example, it has been discovered that the two-loop sunrise 
integral with unequal masses in two space-time dimensions satisfies an ordinary second-order 
linear differential equation [10|. In the equal mass case, such a second-order linear differential 
equation has been known already for long time [6J. The surprise was that the generalisation to 
unequal masses does not increase the order of the differential equation. This can be understood 
in terms of variations of Hodge structures: The order of the differential equation is related to 
the dimension of a certain cohomology group, both numbers being integers. It turns out that the 
variation from the equal mass case to the unequal mass case is smooth and in particular the order 
of the differential equation is a smooth function. It is clear that a smooth function which takes 
integer values has to be constant. In general, Feynman integrals satisfy differential equations of 
Picard-Fuchs type (ordinary linear differential equations with at most regular singularities) [i25l . 
In this paper we solve the differential equation for the unequal mass case, similar to what has 
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been done in the simpler case of equal masses in [|6l. As a side remark, we would like to mention 
that inhomogeneous Picard-Fuchs equations have also been discussed recently in the context of 
mirror symmetry [l26l . 

This paper is organised as follows: In the next section we discuss the relation between (fac- 
torised) differential equations and iterated integrals. The case where the differential operator fac- 
torises into linear factors will lead to multiple poly logarithms. In this paper we go beyond linear 
factors and we are interested in the case, where the differential equation contains an irreducible 
second-order differential operator. To give a simple and concrete example for the approach by 
differential equations we consider in section (|3]) as a warm-up exercise the differential equa- 
tion for the one-loop two-point function and its solution. We then turn to the two-loop sunrise 
integral and briefly recall the second-order differential equation for this integral in section dH). 
Section ([5]) is the main part of this paper. We discuss the solution of the differential equation 
by first determining the solutions of the homogeneous equation, fixing the boundary values and 
by considering one special solution of the inhomogeneous equation. Finally section Q contains 
our conclusions. In an appendix we give useful information on how to transform an elliptic inte- 
gral to the Legendre normal form. We also give a representation of the solution of the two-loop 
sunrise integral in terms of multiple sums. 

2 Motivation 

We denote by Iq a scalar Feynman integral corresponding to a Feynman graph G. A scalar Feyn- 
man integral depends on the external momenta pj only through the Lorentz invariant quantities 

^jk = (Pj+Pk) ■ (1) 

The Feynman integral depends also on the internal masses denoted by m,. We call the set 

S = {sjk.mj} (2) 

the set of kinematical invariants. Let us pick one variable t E S. We will consider the Feynman 
integral as a function of t and write /g(0' the dependence on all other parameters is implicitly 
understood. Recently a systematic method has been developed to find for the Feynman integral a 
Picard-Fuchs equation [[251 . The Picard-Fuchs equation is an ordinary linear differential equation 
of a certain order r. 

Y^p^{t)—Mt) = £^,W/g,.(0, (3) 

7=0 "-^ i 

where /g, denotes simpler Feynman integrals corresponding to graphs, where one propagator of 
the original graph has been contracted. The coefficients Pj{t) of the differential operator and the 
coefficients qi{t) appearing on the right hand side are polynomials in t. (Note that in eq. ([3]) Pj{t) 
denotes a polynomial in t, while in eq. ([T]) pj denotes some external momentum.) 
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In the simplest case the differential operator factorises into linear factors: 

ipA')'^, = (<,,(,)|-fM.))...(«2W^-fft2«))(«,(')^+fc,w). (4) 



7=0 

Let us denote the homogeneous solution of the j-th factor by 

bj{s) 



Mfj{t) = exp ( - / ds^^ I . (5) 

bj{s)/aj{s) is a rational function in s and the integral in eq. © can be expressed in terms of 
rational functions and logarithms. In the case where the Picard-Fuchs operator factorises into 
linear factors as in eq. (jU, the solution of the differential equations is given by iterated integrals 

t ( fi 



Ht) = C,Vi(')+Vl(»/ C.V2(>l) + »l>2(..)/; 



ai(<l)>|/l((l) 1 ^ 'J iJ2('2)V2(<2)"' 

The r integration constants are denoted by Ci, 

In this paper we are interested in cases, where the Picard-Fuchs operators do not factorise 
into linear factors. The simplest cases beyond linear factors have an irreducible second-order 
differential operator. An example for an irreducible second-order differential operator is given 
by 

Let us now consider, what happens if we replace in eq. (H) the j-th factor 

ajit)j^+bj{t) (8) 
by the second-order differential operator 

First of all, a second-order differential operator has two independent homogeneous solutions, 

which we denote by 'y^^j^\t) and 'y\f^p{t). Once these are known, we can replace in eq. © the 
combination 

C,vKO-i)+VKO-i)/^;(^ (10) 
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by 



dt 







ajitj)Wj{tj) 

where Wj{t) denotes the Wronski determinant. The Wronski determinant is defined by 

Wjit) = V/\t)j^Vpit)-^f\t)^^^l\t). (11) 

Cj^^ and Cj^^ are two integration constants. Apart from this substitution the structure of the 
solution in eq. Q as an iterated integral remains unchanged. 

The most important ingredients are certainly the two solutions '^^p{t) and \\f^p (t) of the 
homogeneous second-order differential equation. For a second-order dilferential equation this is 
a non-trivial task, even in the case where Qj{t^, bj{t^ and c j{t^ are polynomials. This is in contrast 
to the linear case, where the generic form in eq. (0 of the homogeneous solution is known and 
an algorithm (partial fraction decomposition) is known to perform the integration in eq. (|5]). Let 
us come back to the example in eq. ([7]). This differential operator has as homogeneous solutions 
the functions 



V^)(0 = ^W, V'^(0 = i^(vT^), (12) 

where K{x) denotes the complete elliptic integral of the first kind 



1 

K{x) = I ■ (13) 

We will see in the sequel that the differential equation for the two-loop sunrise integral has 
homogeneous solutions, which are given by complete elliptic integrals of the first kind times 
algebraic pref actors. 



3 A warm-up exercise 

We would like to make the formal considerations of the previous section a little bit more concrete. 
This is best illustrated by a simple Feynman integral having a first-order differential equation. 
Therefore we consider as a warm-up exercise the one-loop two -point functions in two space-time 
dimensions. The integral is defined by 

B{t) = I (fxh{\-xi-X2)—- 2- (14) 

J [—t)x\X2+x\m\+X2m2 
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The analytic result for this integral is well known: 



Bit) 



Pi{t) 



In 



m 



(nil +m2] 



f + m^-t-^pi{t) 



{nil -m2) 



(15) 



The analytical result is easily obtained by direct integration in eq. (fT4l) . As an exercise for the 
extension towards two-loops we derive this result in the following in an alternative way through 
a differential equation. The differential equation for B{t) reads 



B{t) 



Po{t) 



2 2 

t — nil ~^2- 



A solution of the homogeneous equation is given by 



\|/ {t) = exp 

The boundary value at ? = is given by 

5(0) 



dv 



Pojv) 
Piiv) 



— mj 



Plit) 



(16) 



(17) 



A* 



1 "4 



(18) 



The solution of the inhomogeneous equation can be written in the form 

B{t) = Cx\f{t)+y\f{t)[dv ^^Ifj y 

y pi{v)y^{v) 

With \^(0) = 1 and C = 5(0) we have 

Cv|/(0 = 



(19) 



2 

1 "^2 

:ln^. 

mf 



(20) 

^/Pl{t) niX 

The second term in eq. (fT9l ) is a special solution of the inhomogeneous equation. Direct integra- 
tion leads to 



\|/(f) / dv 



(-V) 

Pi{v)xf{v) 



i\ {m\ -\-m\-t-\- ^/piifj 



(21) 



Plugging the results of eq. (1201) and eq. (|2T)) into eq. (|T9l) we recover the analytical solution in 
eq. (fT5l) . The argument of the logarithm can be re- written as 



m 



j+ml-t+^pi{t) 



m^-\-nv^ — t 



a/ {nil + -t + ^/{ml- m2)^ - 1 
a/ (mi + m2)2 - ? - A/(mi - m2)2 - ? 



(22) 



and we finally obtain 



^J{{ml+m2Y -t){{mi-m2Y -t) \\/{mi+ mjY - f - A/(mi - - 1 

We see that the solution consists of an algebraic prefactor / a/ Piit) and a transcendental func- 
tion ln(...). We note that the algebraic prefactor can already be obtained from the homogeneous 
solution in eq. (|20l) . The point t = (mi + 1712)^ is called the threshold, the point t = (mi — m2)^ is 
called the pseudo-threshold. On the principal branch of the logarithm the integral is finite at the 
pseudo-threshold. This is easily seen from eq. (|23T ): At the pseudo-threshold the argument of the 
logarithm is approaching one, thus cancelling the square-root singularity of the prefactor. 



4 Definition of the two-loop sunrise integral 

The two-loop integral corresponding to the sunrise graph with arbitrary masses is given in 2- 
dimensional Minkowski space by 

S(p^) = .^fi^ilHl L _. (24) 

The corresponding sunrise graph is shown in fig. (H}. In eq. (|24)) the three internal masses are 
denoted by mi, m2 and mj. Without loss of generality we assume that the masses are ordered as 

< mi < m2 < m3. (25) 

The arbitrary scale /j is introduced to make the integral dimensionless. denotes the momentum 
squared. This variable plays an important role in our derivation and it is convenient to introduce 
the notation 

t = p^. (26) 
In terms of Feynman parameters the two-loop integral is given by 

S{t) = p^lj (27) 
a 

with 

^ = —xiX2X3t+[ximi+X2m2+X3m2)'U, 

U = X1X2 +-'C2X3 +X3X1. (28) 

The differential two-form CO is given by 

CO = xidx2 A dx3 +X2dx3 A dxi +X3dxi A dx2- (29) 
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Figure 1 : The two-loop sunrise graph. 

Eq. (|27|) is a projective integral and we can integrate over any surface covering the solid angle of 
the first octant: 

o = {[jci :x2 :jc3] gP^|x,->0,/= 1,2,3}. (30) 

It will be convenient to introduce the following notation: We denote the masses related to the 
pseudo-thresholds by 

^I=mi+m2 — m3, 1^2 = mi —m2 + m^, = —mi +m2 + mj, (31) 

and the mass related to the threshold by 

IJ4 = mi+m2 + m3. (32) 

With the ordering of eq. (|25T ) we have 

A'l < ^^2 < W < A'4- (33) 

We further introduce the monomial symmetric polynomials M^^i^i^ in the variables m^, and 
mj. These are defined by 

^x,...3 = IK)^(^')H)^(^^)H)^('^\ (34) 

a 

where the sum is over all distinct permutations of {Xi , X2, ^13) . A few examples are 

^^100 = ml + ml + ml, 
Mill = mjm2m3, 

AO 47 Zl7 Zl7 zl7 47 

M210 = mim2 + ni2ni2+in2nii+m2mi+m2m2 + mim2. (35) 
We also introduce the abbreviations 

A = 2Miio -M200 = A'iA'2A'3A'4, 

di = —ml + m^ + m^, 62 = — m^ + rn^, 63 = ml + m^ — m^. (36) 
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As shown in [[TOl the two-loop sunrise integral in D = 2 dimensions satisfies a second-order 
differential equation: 



S{t) = ^?p^{t), (37) 



where po(0' P2{t) and P3,{t) are polynomials in t. (We follow here the notation of [[TOl 

and denote the coefficient of jdt^ by P2-j{t).) The polynomial po{t) is given by 

poit) = t{t-p\){t-pl){t-4){t-pl){3t^-2M^oot + A). (38) 

The zeros of the polynomial po{t) correspond to the singular points of the differential equation. 
The regular singular points at ? = 0, at the pseudo-thresholds 1^1,^12 and ^3 as well as the threshold 
/J4 are expected. In addition the roots of the quadratic polynomial 3t^ — 2MiQQt + A define two 
regular singular points. We will comment on these points after eq. (|53T ). The polynomials pi (t) 
and P2{t) read 

Pi{t) = 9?^-32Mioo?^ + (37M2oo + 70Miio)?^-(8M3oo + 56M2io + 144Miii)r3 (39) 

- ( 1 3M400 - 36M3 10 46M220 - 1 24M21 1 ) 

- ( -8M500 + 24M410 - 1 6M320 - 96M3 1 1 + 144M221 ) t 

- (Meoo - 6M510 + I5M420 - 2OM330 + I8M411 - I2M321 - 6M222) , 
P2{t) = 3r5-7Mioo?^ + (2M2oo + 16Miio)?3 + (6M3oo-14M2io)?2 

- (5M400 - 8M310 + 6M220 - 8M211) ? + (M500 - 3M410 + 2M320 + 8M311 - IOM221) . 

In contrast to the polynomial po{t), which is given in factorised form in eq. (|38T ). we do not know 
a factorisation for pi (t) or P2{t) ■ However, we can write both polynomials as a sum of two terms, 
one proportional to the quadratic polynomial 3t^ — 2Mioo? + A and the second proportional to the 
derivative 6t — 2A/ioo of the quadratic polynomial: 

Pi{t) = ^po{t)-2po{t)^\n{3t^-2Mioot+A) , (40) 
dt dt 

P2{t) = -^^PQ{t)j\n {3t^ - 2Mioot + A) 

+^ {3t^ - 2MioQt + A) {3t^ - IMioot^ + (5M200 + 6M1 10) t - M300 -f M210 - 1 8M1 1 1 ) . 
The polynomial P3{t) appearing in the inhomogeneous part is given by 

P3{t) = -2{3t^-2Mmt + Af (41) 

2 2 2 

fft tn fTi 

+2c (f,mi,m2,m3)ln + 2c(?,m2,m3,mi) In— ^ + 2c (?,m3,mi,m2) In— |-, 
/j^ iw^ iw^ 
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with 



c(f,mi,m2,m3) = 

[-2ml + fnl + mf) + {6m\ - ?,m\ - 3m^ - lm\m\ - lm\ml + lAm\ml) 

(f\ f\ f\ AO AO OA OA AO 

—6m" + + Sm^ + 1 lmjm2 + 1 Im^m^ — Sm^m2 — Sm^m^ — 3m2m^ — 3m 



lm\)t 



+ {2m\ 



o Q f\ f\ f\ f\ f\ f\ 

m2 — m3 — 5mjm2 — Sm^m^ + m^mj + m^m^^ + Amjm^^ + 'Xm2m-^ 



i-244|o44 z:44|o422 242 22 4\ /^ox 

+3m^m2 + im^mi^ — om2m^ + 2mjm2m3 — m^m2m^ — mim2m^) . (42) 



We remark that the inhomogeneous second-order differential equation in eq. (|37l) implies that 
S{t) is a solution of a factorised homogeneous third-order differential equation: 



S{t) = 0, (43) 



with 



P3{t) = j/3{t). (44) 



5 The solution of the differential equation 



We will solve the second-order differential equation in eq. (1371) by first considering the corre- 
sponding homogeneous equation. We denote the two independent solutions of the homogeneous 
equation by \\fi{t) and >|/2(0- The solution of the inhomogeneous differential equation is then 
obtained through the method of the variation of the constants. The full solution can be written as 

t 

5(0 = Cm{t)+C2y\f2{t)+p^ ldh—^0^[-y^i{t)y\f2{h^ (45) 



where W{t) denotes the Wronski determinant. In the next sub-sections we will first determine 
the homogeneous solutions (t) and V|/2(0' then fix the constants Ci and C2 from the boundary 
values and finally consider a special solution of the inhomogeneous equation, given by the third 
term in eq. (|45T ). 

The integral S{t) is real in the Euclidean region defined by ? < and > 0. The analytical 
continuation to ? > is provided by Feynman's zO-prescription, where we substitute t ^ t + 
iO. The symbol +/0 denotes an infinitesimal positive imaginary part. We will often encounter 



expressions of the form \/nP- — t. We have 

r-^ f M+ for t < m^, 

[ -iR^ for t > m^. 

We can factorise this expression as 



(46) 



\/m^-t-iO = J (rn- Vt+lO^ (m + y/t + iO^ . (47) 
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Note that 



^—77 . I 77 f ?(v^-«0) for ?<0, 

\ v^ + ?0 for r>0. 



(48) 



In the following we will suppress the /0-part and write simply 



\Jm?--t = Y (m - a/z") (m + a/F) . (49) 

5.1 The homogeneous solutions 

In this section we consider the solution of the homogeneous differential equation 



d 



\\f{t) = 0. (50) 



We will denote the two linear independent solutions of eq. (I50l) by ^i{t) and v/2(0- We recall 
that the Wronski determinant is defined by 

W{t) = x,/i(^)lx,/2(0-x|/2(0^¥i(0 (51) 



and satisfies 



^W{t) = -^Wit). (52) 
dt po{t) 

This equation can easily be integrated and yields 

The integration constant c' will be fixed later. The Wronski determinant vanishes at the regular 
singular points defined by the roots of the equation 3t^ — 2Mioot + A = 0. 

In order to present the solutions of the homogeneous equation we introduce the notation 

(mi-m2)^ im^-^Ttf' [m-i + y/tf (mi+m2)^ 
X\ = 5 , X2 = ^ , JC3 = 5 , X4 = 5 . (54) 

/J^ [A"^ jj"^ 

For t <0 the quantities X2 and JC3 are complex. For t>0 the four quantities are real. If in addition 
to ? > we assume < (mi + m2)^ and that t is small they are ordered as 

<Xi < X2 < X2 < X4. (55) 

The equation 

= {x — xi){x — X2){x — X3){x~X4) (56) 
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defines an elliptic curve. The holomorphic one-form associated to this elliptic curve h clx/y and 
the periods of the elliptic curve are given by 



= 2/- = K 

^/[x^-Xi)[XA-X2) 



X3 






dx 




y 






^3 




^/ 


dx 




y 







X|/2(0 = 21- = =K\ 1. (57) 



(■^3 


~X2] 


1 (JC4 


-Xi) 














|(X4 




{X3 




1 (;C4 


-X2) 



K{k) denotes the complete elliptic integral of the first kind 
The modulus k{t) of (t) is given by 



(58) 



j^.. ^ / (X3 -X2) ix4-xi) _ / 16mim2m3v/F 

(^3--«l)(^4-^2) V (-"l+V^) (-"2 + ^?) (W + V^) (a'4-V^)' 



We recall that the pseudo-thresholds jui, jLi2, P3 and the threshold /J4 have been defined in eq. yl 
and in eq. (|32l) . respectively. The complementary modulus is given by 



^//^^ _ / (x2-JCl)(x4-.y3) _ / (a^I - V^) (^^2 " V^) (a'S " V^) + V^) 

Note that we have the relation 

k{tf + k'{tf = 1. (61) 



We would like to remark that although we started in eq. (p6l ) from an equation which is not 
symmetric with respect to the permutation of the masses mi, m2 and m^, the functions (?) and 
\|/2(f) are symmetric in the masses. We further remark that the elliptic curve defined in eq. (|56l ) 
varies with t, and so do the periods >|/i(?) and \|/2(0- 

The main result of this sub-section is the following: The periods \\fi{t) and \\f2it) defined in 
eq. (l57l) are two independent solutions of the homogeneous differential equation. This is easily 
verified by inserting the functions ^i{t) and V|/2(0 into eq. (|50l) . The solutions have been found 
by generalising the corresponding solutions for the equal mass case in ref. [16] to the unequal 
mass case. The appropriate values for the variables xi, X2, X3 and X4 for the unequal mass case 
have been suggested in [i3i|4i|. The representation in eq. (|57T ) of the homogeneous solutions is not 
unique. To see this in an example we set 



Xi{t) = {x3—X2){x4 — Xi) = 16mim2m3\/t/iJ^, 

X2{t) = {X2 -X\){XA-X3) = (/Ji - \ft) {H2 - Vt) (/Js - Vt) + V^) /A^"^, 

^3 (0 = (-^3 - XI ) (jc4 - X2) = (a/1 + Vt) {m + Vt) (a/3 + Vt) {m - Vt) (62) 



12 



We have 



X1+X2 = X3. (63) 



Then any of the six functions 



V3-^i^f#V ^s-^k(u^ 



V2 = ^i^ A/Tf , \|/4 = ^ i^h/T/ , ^\f6 = ^ KiiJ-^], (64) 



V V ^3 / ' V V ^2 / ' V V ^1 / 



is a solution of the homogeneous differential equation. Of course, an ordinary second-order 
differential equation has only two independent solutions, therefore there must be relations among 
these functions. Indeed one finds ETll 

^3 =\^l -5\|/2, V4 = -^Vl+¥2, ¥5=¥l, ¥6 =¥2, (65) 

where s = sign(Im(Xi/X3)). We have chosen the functions and \|/2 as a basis. With this 
choice both the modulus k and the complementary modulus k' are real and in the interval [0, 1] 
(for < (mi + m2)^ and small positive 0- 

The function is regular at ? = and we find the value 

¥i(0) = (66) 
On the other hand, the function '\\f2 exhibits a logarithmic singularity at the origin. We have 

From the solution in eq. (l57l) we find that the constant c' appearing in the Wronski determinant 
is given by 

c' = -2Tii{f?f. (68) 
In simplifying the expression for the Wronski determinant the Legendre relation 

K{k)E{k')+E{k)K{k')-K{k)K{k') = ^ (69) 

is useful, k and k' are the modulus and the complementary modulus defined in eq. (|59l ) and 
eq. (|60l ). respectively. E{x) denotes the complete elliptic integral of the second kind. (The 
definition will be given in eq. (|9T|).) 
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5.2 The boundary values 

In this section we determine the two constants Ci and C2 from the boundary value at ? = 0. We 
have already seen that the homogeneous solution \\f2 has a logarithmic singularity at t = 0. On 
the other hand the sunrise integral is regular at t = 0. Therefore 

C2 = 0. (70) 

The constant Ci is given by 

\/A 

C, = —,S„, (71) 

where Sq is the sunrise integral at t = 0: 

So = fJ^ 7 9 7 ON — • (72) 

J \x\m\+X2m2 + XT,mi^j U 

We can integrate over any surface covering the solid angle of the first octant and it is convenient 
to do the integration in the plane 

9 9 9 9 

x\mi-\-X2m2-\-XT,mT^ = M , (73) 

where is an arbitrary constant. The change of variables 

222 
, mf , mi .mi 

^i = ^^i, ^2 = ^^2, -^3 = ^^3 (74) 
transforms the integral into 

So = — r ^ . (75) 

J m + m2X3X\ + mi^x\X2 

This is nothing else than the Feynman parameter integral for the one-loop three-point function in 
four dimensions with massless internal lines and three external masses fl28l . With P1+P2 + P3 — 
Oand 

T{4,Pi,P2,P3,P ) = P / — 2 ; —7 

J ni\.k^)^-^k-pi)^)[-{k-p,-p2) 

^ f ^ (76) 

I i-pi)^2X3 + {-pl)x3Xi + {~pl)xiX2 

we have 

So = T{4,-mi-ml-ml^^). (77) 
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In a neighbourhood of the equal mass point we have A > and the integral Sq is expressed in the 
region A > by n29ti3Tl 

5o = ^ [CI2 (ai) +CI2 (as) +CI2 (as)] . (78) 
VA 

The Clausen function Cl2{x) is given in terms of dilogarithms by 

CI2W = i[Li2(.")-Li2 (79) 
The arguments a, of the Clausen functions are defined by 

a,- = 2arctan^^j, /G {1,2,3}. (80) 

We recall that A and the quantities 5, have been defined in eq. (|36l ). Let us summarise what we 
obtained so far: The homogeneous solution with the correct boundary value at ? = is given by 

CiVl(0+C2V|/2(0= (81) 

^^^^^^^ [CI2 (ai) +CI2 (a2) +CI2 (as)] . 

{/Ui + y/t) {lU2 + Vt) (a^3 + Vt) (a^4 - Vt) 

5.3 The inhomogeneous solutions 

In this section we consider a special solution of the inhomogeneous equation. Variation of the 
constants leads to 

f 

SspecialW = dti ffllll [-yi(0x|/2(?l) + X|/2(0¥l (^l)] • (82) 

i Po{ti)W{ti) 

Note that 

2 P3{t) 1 1 1 



(83) 



Poit)W{t) infi^ m'^ {3t^-2Mioot + Af 

c mi , m2, ms) In ^ + c (f , m2, m3, mi) In — ^ + c (f, m3, mi , m2j In 

jJ.'^ /J^ 

The function c mi , m2, ms) is a polynomial of degree 3 in the variable t. Let us denote the roots 
of the quadratic equation 



3t^-2MiQQt + A = (84) 



by f+ and . We have 



t± = ^ (mioo ± 2 v/M2oo - Mho) • 



(85) 
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Partial fraction decomposition leads to 

.2 P3(t) 1 



po{t)W{t) 



^ ^ a2 (?+,mi,m2,m3) _^ a2 [t^ , mi , m2, ma) 



ai (?+,mi,m2,m3) ai (f_,mi,m2,m3) 



(86) 



with 



ai{t,mi,m2,m3) = Ci{t,mi,m2,m3)\n^ +Ci{t,m2,m3,mi)\n^ +Ci{t,m3,mi,m2)m^ 

//^ //^ //^ 



and 



ci (?,mi,m2,m3) = - [—2mi + m2 + + mim2 + mim^ — 2m2m^ + t (2mi — m2 — ma)] , 
C2{t,mi,m2,m3) 



[A(2m| — m2 — m3 — 2mf m2 — 2mf m3 + 4m2m3) 



9(M2oo-Miio) 

+t{2mi — m2 — m3 +9mim2 + 9mim7, — 6mjm2 — 6mjm3 — 3m2m3 — 3m2mj) 



(87) 



We therefore have 



^special (0 = ^^^1 [¥2 (OVl (?l) "¥1 (0^2 (^l)] 



(88) 



a2{t+,mi,m2,m3) a2{t-,mi,m2,m3) ai (r+,mi,OT2, ^3) ai {t^,mi,m2,m3) 



(^i-^-r 



(?+-?_)(ri-f+) (r+-r_)(ri-r_) 



We expect the poles at r± to be non-physical and spurious. They can be eliminated as follows: 
We recall that the functions (t) and \|/2(0 have been defined as 



¥1 (0 



--K{k{t)), X|/2(0 



4/ 



Let us define two associated functions (|)i (t) and (|)2(0 by 



K{k'{t)) 



(89) 



^1 (0 = -^j=^ [K{k{t))-E{km, M) = it)) 
where £'(jc) denotes the complete elliptic integral of the second kind: 



(90) 



E{x) = I 



dt- 



(91) 
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The four functions ^i{t), ^2(0' ^^'^ ^2(0 the entries of the period matrix P of the 
elliptic curve in eq. (1561 ): 



(92) 



There is an algebraic relation between the four functions v|ri, v|r2, and <^2- From the Legendre 
relation it follows that 

XKl (0(^2(0- ¥2(0(^1(0 = TTT-T- (93) 

In order to eliminate the double poles l/(?i — ?±)^ we start from the identity 

' (94) 



(?i -?±)^ dti \ti-t± t± 

and use integration-by-parts. The constant term on the right-hand side of eq. (|94l) ensures that 
boundary terms vanish. The derivatives of \|/2, and (|)2 are given by 

d I d ^ Id X2 

j^^i = -2V'-^l"^2+2C^,^ln-, 

d ^ I d X2 Id X2 

We recall from the definition of Xi{t) in eq. (|62|) that these functions factorise in the variable y^. 
As a consequence, 

^lnX,(0 (96) 
dt 

has a partial fraction decomposition in \/t. We denote 

r[i{h) = V2(0¥i(2'i)-¥i(0¥2(?i), 

ri2(?i) = ^2{t)(^i{h)-\\ri{t)(^2{h). (97) 

Using eq. (|94l ) for integration-by-parts in combination with eq. (|95] ) one finds that this procedure 
does not only eliminate the double poles l/(?i — ?±)^, but also that the coefficients of the single 
poles l/(?i — 1±) are zero. We obtain 

1 /" , r / N Z?! (mi,m2,m3)?i — Z?o('Wi,m2,m3) , , , , .A 

(98) 

with 

222 
tyi tyi fn 

bi {nil, 1712, m^) = di {mi,m2,m^) In^ + di (m2,m3,mi) In^ +di {m^,mi,m2) In^ (99) 

fj^ jj^ jj-^ 
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and 

d\{m\,m2,m2,) — 2m\ — m^ — m^, 

A A A 00 00 00 

do{mi,m2,m3) = Im^ — m2 — m3 — mjm2 — mjm3 + 2m2m3. (100) 

Note that in the equal mass case m\=m2 = the terms involving r\2{h) — 'ni(^i) vanish, due 
to the fact that 

di{m,m,m) = 0. (101) 

5.4 The full result 

For the convenience of the reader we put here all pieces together and recall the relevant def- 
initions. The complete result for the two-loop sunrise integral with arbitrary masses in two 
space-time dimensions is given by 

S{t) = - [Ch (ai ) + CI2 (a2) + CI2 (as)] x|/i (t) 

1 /„ N bi{mi,m2,m3)ti-bo{mi,m2,m3) . .\ 
-2 J dH |Tii {H) Ji^i;) ["12 {h) -Til {h)] I , (102) 



inn 





where 



rii(?i) = V2(0¥i(^i)-¥i(0¥2(?i), 

V[2{t\) = \/2(04'i(^i)-¥i(0^2(/'i). (103) 
The functions \|/2, ^1 and ^2 are given by 

^1 = ^^)) - ^ ^^))] ' ^2 (t) = -^=E {k' (0) . (104) 

The modulus k{t) and the complementary modulus }<f{t) are defined by 



ao5) 



The functions Xi(?), X2{t) andX3(?) are defined by 

Xi{t) = \6m\m2m'i\/t / ^ , 



X2 (t) = {ni - Vt) {m2 - Vt) {n3 - Vt) {^4 + Vt) //, 

X3{t) = {iJi + y/i){iJ2+Vi){n3 + Vi){m-Vi)/ij\ (106) 
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i"!, jU2, jU3 and //4 are defined by 



= mi +m2 — m3, = mi — m2 +m3, = — mi +m2 +m3, = mi +m2 +m3. (107) 
The angles ai, a2 and are defined by 

a; = 2arctan j ' (10^) 

with 

ooo ooo ooo 

A = jui^2j"3j"4, 5i = — mi +m2 + m3, 52 = mi— m2+m3, 53=mi+m2 — m3. (109) 
Finally, the functions &i(mi,m2,m3) and bo{mi,m2,m2) are given by 

2 2 2 

TYl fYt f7% 

bi (mi , m2, m3) = J( (mi , m2, m3) In ^ + di (m2, ms, mi ) In ^ + J, (m3 , mi , m2) In ^ , (110) 

//^ //^ //^ 

with 

Jl (mi,m2,m3) = 2mi— m2 — m3, 

dQ{mi,m2,m2) = 2mi—m2 — m^ — m^m2 — m^m^ + 2m2rn^. (HI) 

6 Conclusions 

In this paper we discussed the analytic solution of the two-loop sunrise integral with unequal 
non-zero masses in two space-time dimensions. We obtained the solution by solving an inhomo- 
geneous second-order linear differential equation. The homogeneous solutions involve complete 
elliptic integrals of the first kind. The two-loop sunrise integral with non-zero masses is the 
simplest Feynman integral, which cannot be expressed in terms of multiple polylogarithms. The 
complications are related to an irreducible second-order differential operator. We expect the 
methods discussed in this paper to be useful for other computations as well. 
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A Transformation to Legendre normal form 

In this appendix we discuss how to transform the general quartic elliptic curve 

= (x — xi) (x — X2) (x — X2) (x — X4) (112) 
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to the cubic Legendre normal form 

w2 = (113) 
Without loss of generality we can assume that 

Xi<X2<X3<X4. (114) 

The mapping 

{X3 -X4){x-Xi) XI {X4 - X3) + X4 ix3 - Xi) Z 



Z = 



{X3-Xi){x-X4)' X4-X3 + {X3-Xi)z 

dx (X3 — ;ici) (;C4 — X3) (;C4 — Xi' 



dz [X4-X3 + {X3-Xi)z\^ 

sends the points x\,X3, X4 to the points 0, 1 and infinity. The point X2 is mapped to 



We thus have 

dx \ dz \ dz 



y ^J{x3-Xi){x4-X2) W a/(.«3--^i)(-«4--^2) a/z(1 - z) {z-X) 

and with the notation of eq. (|60l) k' = a/X. With this mapping we obtain for the integrals 



X3 1 
dx 1 f dz 



y \/{x3-X\) {X4-X2) J W a/(x3-Xi) {X4-X2) 

1 

dx 1 f dz 2i 



I 



X4 



y y^{x3-xi) {X4-X2) J w ^j[^^zr^^^y(^^zr^ 



(115) 



(x2-.yi)(jC4-JC3) 

k = -. (116) 

(;c3-;ci)(;c4-;c2) 



(117) 



, Kjl). (118) 

J w ^/(x-^ — xi) (xa — x^) 



B Representation in terms of Lauricella functions 

The main focus of this paper is to represent the two-loop sunrise integral as an iterated integral 
in the spirit of section |2l Besides iterated integrals a representation in terms of nested sums can 
also be useful. In this appendix we briefly review a nested sum representation for the two-loop 
sunrise integral. In D space-time dimensions this representation is given in terms of Lauricella 
functions. Unfortunately this representation does not allow a straightforward specialisation to 
two space-time dimensions. This is due to spurious singularities which appear in intermediate 
results. We will discuss this issue in the following. We start with the two-loop sunrise integral in 
D space-time dimensions. 

/ 2^3-^ fd^ki d^k2 1 
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The first step is to use the Mellin-Bames technique to write the propagators as 

■ + 



1 



1 r im^V 

— I dor(-(5)r(o+i) ^ i,, . 



(120) 



The integration in eq. (11201 ) is over a contour along the imaginary axis separating the poles of 
r(— o) and r(o+ 1). The integration over the loop momenta ki and ^2 corresponds then to an 
integration for a massless sunrise integral and can be done easily. Setting 



m, 



mi, 



Xl 



X2 



i-ty i-t) 

one arrives at the Mellin-Bames representation 



X3 



(121) 



5(0 



tx""-^ 1 

— I JOi JO2 J03 4'x2'4T(-Oi)r(-02)r(-03) 

(271/) 



r{(5i+02 + 03 + 3-D) 



r( -oi + ^-i )r(-<52 + ?-i )r(-<'3 + ?-i < , , 

2 ; V 2 ; V 2 ;r(-oi-02-03 + §D-3) 

The next step is then to close the contours and to take residues. For \xi\ < 1 (/ = 1,2,3) we can 
close all contours to the right and we obtain a symmetric result in xi, X2 and X3. The result can 
be expressed in terms of the third Lauricella function in three variables, defined by 



cx) 00 00 



Fc{ai,a2;bi,b2,b3;xi,X2,X3) = H H 

mi —Q in2~0 m2~0 



(^l)mi- 



/ \ nil '"2 "'3 

m2+m3 ('^2jmi+m2+m3 -"-J -^2 '^3 



1112 



mi ! m2 ! m3 ! 



(122) 



where (a)„ = r{a + n) /r(a) denotes the Pochhammer symbol. We obtain 

_ D-3 

r(3-D)r(^-i)^ ^ 3 D D D 

x< r(4-3) -^c(3-A4--Z);2--,2--,2--;-Jci, -..2,-^3) 



r(2-f)r(i-f)r(f-iy 



r(D-2) 



D ^ D 



Fc{3 - D, 2 - -; -, 2 - -, 2 - -; -xi, -X2, -X3).x:i2 



D 



+Fc(3 - D, 2 - -; z - -, -, z - -; -jci , -JC2, -JC3j^f 



+Fc(3 -D, 2 - -; 2 - -, 2 - -, -; -jci , -.)C2, -X3)x| 



-1 



-1 



D D D 



D 



^c(l,2 2 --;-xi,-X2,-JC3) (■«i-«2 



+Fc{l,2-—;—,2-—,—;-xi,-X2,-X3){xiX3)2 

,^ ^ D ^ D D D , , , D 1" 

+Fc(l,2--;2--,-,-;-xi,-.x:2,-JC3) (^2^3)^ 



(123) 
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This expresses the sunrise integral in D space-time dimensions as a linear combination of seven 
Lauricella functions. This result can be found already in [|2l|, however eq. (24) of [|2l has a few 
typos. For the convenience of the reader we give here the correct result. 

Eq. (|123l) expresses the sunrise integral in D space-time dimensions as a linear combination 
of seven Lauricella functions. We are interested in the result in two space-time dimensions. In 
two space-time dimensions, eq. (11231) is not yet particular useful, because the individual terms 
diverge in the D — )■ 2 limit due to the prefactors Y{D /l—l) and r( 1 — D/2) . If we set D = 2 — 2£, 
then each of the seven terms has a Laurent series in £ with poles up to l/£^. The poles are 
spurious singularities and cancel in the sum. We obtain a manifest finite result for two-space 
time dimensions by expanding in £. In order to present the result, we introduce the Euler-Zagier 
sums (321 

Zi(n) = £i Zii(n) = £lzi(j-1). (124) 
Then the result is given by 

2 oo OO OO / • |\2 

^ > jl=0j2=0j3=0 \J^-J2-J^-/ 
X {12Zii (Jl23) +6Zi (ji23)Zi (Jl23) -8Zi (Ji23) [Zi (ji) +Zi (72) +Zi (73)] 

+4[Zi (7i)Zi (j2) +Zi (j2)Zi (j3) +Zi (j3)Zi (ji)] +2 [2Zi (7123) -Zi (72) -Zi (73)]lnxi 

+2[2Zi Um)-Zi {h)-Zi (ji)] ln.^2 + 2 [2Zi (7123) -Zi {ji) -Zi (72)] lnx3 

-(-Inxi lnx2 + Inx2lnx3 + lnx3lnxi} . (125) 

Here we used the notation J123 = ji + J2 + J3- This expresses the two-loop sunrise integral in 
two space-time dimensions as a five-fold sum. 
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